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^^ Abstract. These proceedings consist of a brief overview of the current understanding 

^ of collective behavior in relativistic heavy-ion collisions. In particular, recent progress 

^^ in understanding the implications of event-by-event fluctuations have solved important 

'^ puzzles in existing data — the "ridge" and "shoulder" phenomena of long-range two- 

,—1 particle correlations — and have created an exciting opportunity to tightly constrain 

1 -p^ theoretical models with many new observables. 

I — I 1. Introduction 

^ It is well established that the medium created in relativistic heavy-ion collision is 

^ characterized by strong collective behavior. In particular the large "elliptic flow" 

v/~j anisotropy in multiparticle correlations has been interpreted as evidence of a low- 

viscosity, nearly-perfect fluid being created in these experiments. 
f — . 
(^ Even so, it was only recently that we have begun to appreciate the full impact of 

'""' this collective evolution on these correlation measurements, due to the non-trivial effect 

L| of event-by-event fluctuations. In these proceedings I review this recent progress and 

speculate on future prospects for new flow measurements to deepen our understanding 
of all stages of the evolution of heavy-ion collisions. 
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2. Two-particle correlations 

Correlations between particles detected in heavy-ion collisions have been a powerful tool 
in understanding the properties and dynamics of the collision system. The simplest and 
most ubiquitous is a two-particle correlation, where the distribution of pairs of particles 
from the same event are calculated and averaged over events. Since the azimuthal 
orientation of each collision is uncontrolled, one can only measure rotationally symmetric 
quantities. Therefore, one must fix the relative azimuthal angle of the particle pair 
A0 = (f)^^^ — (j)^"^^ and measure the average dNpi^ij-s/dAcj) (the pseudorapidity dependence 
is not so constrained, but is also often presented as a function of Arj). Two recent 
examples are shown in Fig. [T] 
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Figure 1. Preliminary two-particle correlation data for charged hadrons in Pb-Pb 
collisions from the CMS (left) [1] and ATLAS (right) [2] collaborations as presented 
at this Quark Matter 2011 conference. The CMS analysis restricts one particle of the 
pair to have transverse momentum 2 < pt < 4 GeV and the other 4 < pt < 6 GeV 
while ATLAS uses all particles with 2 < pt < 3 GeV. 



At small relative pseudorapidity are features that are familiar from smaller collision 
systems like proton-proton — e.g., the positive correlations near A0 = At] = due to 
Bose-Einstein correlations and jets. However, a significant portion of the measured 
correlation is unique to heavy-ion collisions and has been interpreted as being generated 
from strong collective motion of the system, or "flow" — in particular the large second 
Fourier component in A0 that is present even for particles separated by a very large 
relative pseudorapidity Arj. 

The standard picture of this flow correlation is as follows. One imagines that 
particles in a given collision event are emitted at the end of collective evolution of the 
system according to a probability distribution (which depends on the initial geometry 
and subsequent collective expansion as described below). The azimuthal dependence of 
this final particle distribution can be generally written as a Fourier series with respect 
to the azimuthal angle (j) of the momentum of the outgoing particle: 

C) J AT CO OO 

^^^^' =l + ^2^;„cosn(0-v[/„)= ^ t;„e»*"e— <^. (1) 



N 



n=l 



Here Vn are the fiow coefficients, while \&„ is the reference angle for each harmonic 
n that is defined as the phase of the complex Fourier coefficient {e^'^'^) = f„e*"*", or 
equivalently the angle defined by (sinn(0 — \E'„)) = 0. Note that, when defined this way, 
both \l/„ and f„ can in principle depend on transverse momentum and pseudorapidity. 
If there is strong collective behavior, each particle is emitted independently and 
the pair distribution is given simply by the product of single-particle distributions. For 
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example, a single Fourier component of the pair distribution is given by 

where the small inner brackets represent an average over pairs in a single event and the 
large outer brackets indicate an average over events in a centrality class. The entire 
azimuthal dependence can then be written as 

^{^) "^> 1^1 2 (.<>.?) co.„(A,). (3, 

This form is typically called a "flow" correlation, while any intrinsic pair correlations 
in a given event break this single-event factorization, and such an extra contribution 
to the observable can be termed "non-flow". It should be noted that, in addition to 
independent emission, this commonly-written form of Eq. ([s]) also contains the implicit 
assumption that the orientation angle \&„ is a global quantity that depends little on 
transverse momentum or pseudorapidity (or more precisely, ^^^^ ~ \E'^^^ for all pairs in 
the analysis), as well as the assumption that multiplicity fluctuations in the centrality 
bin are small enough that one can take the iVpairs factor out of the average. 

Until recently it was assumed that, due to the symmetry of colliding identical 
spherical nuclei, the following are true: \E'„ ^ when defined with respect to the 
impact parameter, the flow coefficients for ra=odd are odd in (pseudo) rapidity, and 
therefore only even terms n = 2,4 are important near mid-rapidity. Thus, the sizeable 
second Fourier harmonic, largely independent of relative pseudorapidity, that is seen in 
these two-particle correlations has been interpreted as coming from the "elliptic flow" 
response to an initial almond-shaped collision overlap region, and indeed serves as a 
measurement of the coefficient V2- One can see this by eye in Fig. [T] as the two large 
elongated bumps centered at Acj) = and A0 = n. Upon closer inspection, however, it 
becomes apparent that the peaks are not identical. Thus, the A0 dependence can not 
be entirely described by only a second Fourier component, even when gaps are imposed 
in the relative pseudorapidity of the particle pairs to suppress most known non-flow 
correlations. 

The remaining long-range (and until recently, presumed to be non-flow) 
contribution has been a topic of signiflcant study. To do this, one must choose a model 
scheme to separate the flow from non-flow. A common example is the Zero-Yield-At- 
Minimum (ZYAM) prescription |3], which posits a small non-flow correlation coming 
from a set of particles that have zero production at one or more values of A0. Next, 
one must also posit the value of each flow coefficient f„, though in most cases only a 
non-zero V2 is used. After this flow-subtraction procedure, there typically remains a 
narrow "ridge" feature at A0 = 0, as well as a broad feature on the away side, which 
often displays a "shoulder" or "cone" structure with a dip at A0 = vr. Similar features 
can also now be seen in extremely central collisions without subtraction [21 II]- 

Any event-by-event fluctuation, however, will break the apparent symmetry of the 
collision system. Flow fluctuations have been known to be important for several years, 
but initially only the effect on elliptic flow V2 was studied [S] . A crucial step forward was 
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the realization that event-by-event fluctuations imply that other coefficients, especially 
V3, should be non- negligible [6]. Thus, the previous flow-subtraction procedures did not 
actually remove the flow correlation, and in fact, flow correlations alone can potentially 
generate the entire observed two-particle correlation at large At] [7] . There is a growing 
consensus that this is indeed the case, but to test this one must flrst specify all the 
properties we expect such flow correlations to have, and then search the data for evidence 
of a non-flow contribution. 

3. Hydrodynamics 

The state-of-the art description of the bulk evolution of heavy-ion collisions is the 
hydrodynamic framework [8]. Such models characterize the collision medium as an 
approximately thermalized relativistic fluid with a few macroscopic properties such 
as viscosity and equilibrium equation of state, subject to initial conditions and a 
prescription for particles to "freeze-out" from the fluid at the end of the evolution. 
The generic picture is that the spatial anisotropy in the density of the system at the 
onset of collective evolution is converted into a momentum anisotropy that manifests 
itself in correlations among the flnal particles. Event-by-event fluctuations in the flnal 
distribution ([I]), then, ultimately stem from fluctuations in the initial geometry, which 
are evolved forward in time in a hydrodynamic framework. 

The early-time, pre-equilibrium dynamics are not yet understood in detail. 
Typically a few simple models are used to probe the range of results that can reasonably 
be expected. However, there are generic features that are expected. For example, most 
models predict flux-tubes or string-like structures to appear at early times that are 
extended longitudinally. These structures ultimately cause very long-range correlations, 
like those seen in data (see Fig. [I]). The size of event-by-event fluctuations, and thus 
harmonics like ^3 that are generated by fluctuations, should also vary with centrality 
in a predictable way due to the change in size of the system. In addition, the emission 
of approximately thermalized matter at the end of hydrodynamic evolution generically 
results in a mass ordering of each flow coefficient, like that seen experimentally for V2, 
with smaller values at a fixed pt for heavier particles. 

The hydrodynamic response to the initial geometry has been extensively studied 
for smooth and symmetric initial conditions, and more recently with event-by- 
event fluctuating initial conditions [9l [TOl [HI [121 [131 IH] ^^ addition to transport 
calculations fTEl [TB] . Thus, much is already known. However, work is ongoing to 
more precisely characterize this (non-linear) hydrodynamic response in order to more 
precisely extract properties of the collision system. To that end, I describe a systematic 
approach, inspired by Ref. [IT], which more clearly explains the reason behind simple 
approximate hydrodynamic response relations that have been found such as V2 oc 82, 
and suggests how to systematically implement corrections. 

Take the 2D Fourier transform of the initial transverse density p(x) and expand in 
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harmonics and in powers of fc = vk^: 



with 



and 



p(0)- /rf2xp(x) A^o-' ' ' ^^ 



777 ' 

{rV^}, (5) 



2ml m+n \\( m-n \\ 



^ J rf^xp(x) . . ■ 

where now is a spatial coordinate. Note that the moments em,n are only non-zero for 
m > |n| and {m — n) = even. Equivalently, and perhaps more naturally, one can do 
the same decomposition with the log of the Fourier transform to obtain the cumulants 
Wn^m [IZj; but the higher order expressions are more complicated to write. 

The final particle distribution at a given momentum can then be regarded as simply 
a function of (an infinite number of) moments of the initial density distribution. This 
is useful for two reasons. First, hydrodynamics is a long-wavelength description of a 
system that is only valid if short-distance dynamics is not important. Thus, we expect 
the hydrodynamic response to be insensitive to short-distance features represented by 
moments with high powers of fc — i.e., viscous effects dampen modes with larger m. We 
can then expect to be able to truncate this series and write a hydrodynamic response 
function of a finite number of moments, and if the anisotropy is not large compared to the 
system size, we can also write it as a Taylor series with a finite number of terms. Second, 
it is a nice way of clearly specifying the symmetry constraints — each flow coefficient Vn 
can only depend on aspects of the initial density with the correct symmetry properties. 

Take as an example elliptic flow. The second complex Fourier coefficient (e*^*^) = 
y^gi2*2 (.g^j^ only depend on combinations of moments that are dimensionless and have the 
correct symmetries — e.g., a rotation of the system — )■ 0-1-5 results in a multiplicative 
factor e*^'', and if there is a reflection symmetry about some plane in the initial state, 
it must also be present. Choosing a coordinate system such that {x} = {y} = 0, the 
lowest order moments are £:2,o(oc {r'^}) and £2,2(0!: {r^e^*'^}). Thus, to lowest order in m, 
V2 can only depend on the dimensionless ratio 

where epp and ^pp are the standard participant eccentricity and participant plane, 
respectively. Since typically epp ^ 1, we can write the function as a power series in 
epp. Thus, to flrst order we have 

r»,2 2i</.| 
(e*2<^) = V2e'^''' = C eppe^**^^ = -C ^,0. % (8) 

|^2| 

for some (real) coefficient C that encodes information about, e.g., viscosity, equation 
of state, etc. This is the familiar relationship V2 oc £2, which has long been known for 
hydrodynamic simulations with smooth initial conditions, but also works to a reasonable 
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Figure 2. Schematic diagram illustrating the simultaneous orientation of directed 
(^i), elliptic (^2), and triangular flow {"^3) in relation to the initial distribution of 
participant nucleons in a single event from a Glauber Monte Carlo [T5]. 



approximation in lumpy, event-by-event hydrodynamic calculations [9l [TOl [TT] . In this 
context, one can see that the expression is just the first term in a controlled expansion, 
with corrections coming from terms higher order in m (£4,2 oc {r^e^*'^}), or in the Taylor 
series (spp). 

As an aside, it should by now be clear that f 2 can not depend on a term that is 
linear in £3^3, as proposed in Ref. [10], because it does not have the correct symmetries. 
It would have to depend on combinations like £3,353 1 or el^el2\, etc. 

Similar approximate proportionality relations have been found to reasonably well 
describe the results for ^3 [13 E] and Vi [12], while v^ and v^ are more complicated [TT] . 
In retrospect, this is unsurprising since the possible V4 terms £4,4 and £22 ^^^ typically 
of the same size, with the former being more important in central collisions and the 
latter more important in peripheral collisions, in agreement with results from Ref. [TT] . 
Explicitly this could read something like: 



{e'^'t') = v^e 






(9) 



A similar statement can be made about the dependence of v^ on £5^5 and £3,3£2,2- 

The hydrodynamic response has been confirmed to significantly damp higher 
harmonics [18], in agreement with data [2j. Thus, once the hydrodynamic response is 
mapped out for the first ~6 flow harmonics to the order desired, for each centrality and 
each set of parameter values, all useful information about the hydrodynamic model is 
known. This makes it clear exactly what properties of the initial geometry are important, 
and allows one to quickly calculate correlations arising from an arbitrary set of initial 
conditions. 

4. Flow vs. data 



Now that we have a picture of flow, one can look in detail at the long-range two-particle 
correlation data to see whether they quantitatively agree with this picture, or if one 
should instead conclude that other correlations are likely to be present. 
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Already mentioned are a number of generic qualitative features that are in 
agreement with this flow picture. The observed correlation extends to large values 
of A77, with all Fourier components having little dependence on relative pseudorapidity, 
in agreement with flow resulting from the longitudinally-extended structure of most 
models of the initial state. (Note that the flow picture does not predict or require 
exactly zero dependence on Arj). In addition is the lack of non-zero high Fourier 
components, consistent with the observation of strong damping of higher modes in 
hydrodynamics (especially with a non-zero viscosity). Finally, the mass ordering 
predicted by hydrodynamics and seen for V2 is also present in the other Fourier 
harmonics |1]. 

More convincing, of course, is the quantitative agreement with predictions. The 
centrality dependence and size of the second and third harmonics agree with that 
predicted for V2 and ^3 |6l [18]. The pr dependence of all the harmonics as well as 
how they vary with the orientation of the particle pair with respect to the measured 
event plane ^2 is also exactly as expected from flow correlations |7j , with the exception 
of the first harmonic, which non-trivially displays just the prdependence expected from 
a combination of flow plus the well-known momentum-conservation correlation that 
naturally arises when a flnite number of flnal-state particles share a flxed total transverse 
momentum and only contributes to the flrst Fourier harmonic [20] (note that having 
zero transverse momentum does not mean fi = 0, only that [pt ■ Vi) integrates to zero). 

Finally, this Quark Matter conference has seen an even more stringent test — the 
factorization of each Fourier harmonic into a function of the transverse momentum of 
the flrst particle times the pt of the second; i.e., an ansatz something like 

{V^n^vl^^) = V^:^miPt) X V^:^2}{p,). (10) 

Note that this form is even more restrictive than the independent emission assumption of 
a flow correlation, which implies only single-event factorization. For this to be true after 
averaging over events requires that the event-by-event fluctuation can be represented by 
a global factor for each harmonic in each event — e.g., the pt dependence is similar in 
every event and only scales by a constant factor event-by-event. Although this feature 
is not generically required in the flow picture, and hasn't yet been thoroughly tested 
in hydrodynamic calculations, it is a reasonable expectation and a very non-trivial test 
of the existence of non-flow correlations like an away-side jet correlation. Remarkably, 
this factorization has been found to hold quite well [21] for the bulk of particles and for 
every harmonic except the flrst, which is contaminated by the momentum-conservation 
correlation mentioned above. One only sees the factorization begin to break down for 
particles above ~4-5 GeV, and for very peripheral collisions, as expected. 

Hence the growing consensus that, aside from the exceptions listed, flow correlations 
from a hydrodynamic response of fluctuating initial geometries are responsible for the 
entire two-particle correlation structure at large A77. 

A comment should be made here about the so-called "event plane" flow analysis, 
which is also often used to extract flow coefficients f„. Exactly what is measured 
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depends on a quantity called the "event plane resolution", which is essentially the 
square root of the number of particles used in the analysis times the magnitude of 
the flow coefficient itself VnvN [22]. In most cases, the analysis is sensitive only to two- 
particle correlations. If the event plane resolution parameter is very large (as for some 
of the new measurements of elhptic flow at the LHC), the result begins to be sensitive to 
higher particle correlations. The difference between the result of an event-plane analysis 
and a two-particle correlation analysis in such a case (like the difference between two- 
particle and multiparticle analyses) thus depends on non-flow correlations that may be 
present as well as fluctuations, but only in a particular combination that cannot be 
disentangled [22]. One must therefore use caution when comparing w„ from an event 
plane flow measurement to the same quantity extracted from a two-particle correlation, 
in an effort to make a statement about non-flow correlations. If the two values agree, it 
does not indicate an absence of non-flow correlations, only that the event plane resolution 
is not large. If the values do not agree, it could be due entirely to flow fluctuations, 
and similarly this does not indicate the existence of non-flow correlations. In neither 
case does it address the question of this section, whether long-range correlations are due 
only to flow — in other words, taking Vn from a dedicated flow analysis, plugging them 
into Eq. (jsl), and comparing to the total two-particle correlation does not necessarily 
test the presence of non-flow. Likewise, this means that there exists no independent 
measurement of the flow coefficients that one can use in subtraction schemes such as 
ZYAM, and this is why assumptions are always necessary. A general and rigorous way 
to separate flow from non-flow correlations does not exist. 

5. New flow observables 

With the knowledge that long-range two-particle correlations are likely dominated by 
flow correlations comes the opportunity for precise measurement of many independent 
flow observables with little non-flow contamination. The new measurements of v^, W4, 
f5, and vq presented at Quark Matter 2011 [211 [211 121 1251 126] will undoubtedly guide 
theoretical progress for some time, and a dedicated analysis of a rapidity-even component 
of vi [201 [T7] will exhaust the possibilities of strictly two-particle correlations. Already 
from simultaneous measurements of V2 and ^3, one can rule out models for the initial 
conditions (e.g., a particular MC-KLN implementation seems to have insufficient relative 
eccentricity fluctuations) [231 [18] . 

However, an even larger set of observables opens up with correlations between three 
or more particles. Some are familiar, like the fourth cumulant elliptic flow V2{4}, which 
can be generalized to all harmonics, like the recently measured f3{4} [21]. Because of 
the existence of fluctuations, this measures independent information about the event- 
by-event distribution of Vn- However, many other observables are possible, giving 
information about, e.g., correlations between the various orientation angles \E'n [27], 
and providing a wealth of independent information. For example, the comparison of 
^2(4} compared to a two-particle correlation measurement f2{2} had already indicated 
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the above conclusion about the size of fluctuation in the same MC-KLN model [27j . 

With this information from many independent measurements, it should be possible 
not only to more precisely characterized properties of the thermalized medium, such as 
viscosity, but to significantly constrain the nature of fluctuations and the properties of 
the pre-equilibrium system in general. 

6. Summary 

In these proceedings, I have reviewed the emerging realization of the importance of event- 
by-event flow fluctuations in heavy ion collisions. With this knowledge, long-range two- 
particle correlations (e.g., "Mach cone" and "ridge" structures) are naturally explained 
as a result of a collective response to an initial density distribution that fluctuates 
event-by-event. This also implies many new flow observables that can be precisely 
measured and used to significantly restrict theoretical models — not only other Fourier 
components of two-particle correlations like triangular flow, but a wealth of correlations 
between three or more particles — and which will provide exciting opportunities to 
improve our knowledge of the entire collision evolution. 
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